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ABSTRACT 



Context. A large amount of magnetized plasma is frequently ejected from the Sun as coronal mass ejections (CMEs). Some of these 
ejections are detected in the solar wind as magnetic clouds (MCs) that have flux rope signatures. 

Aims. Magnetic clouds are structures that typically expand in the inner heliosphere. We derive the expansion properties of MCs in the 
outer heliosphere from one to five astronomical units to compare them with those in the inner heliosphere. 

Methods. We analyze MCs observed by the Ulysses spacecraft using in situ magnetic field and plasma measurements. The MC 
boundaries are defined in the MC frame after defining the MC axis with a minimum variance method applied only to the flux rope 
structure. As in the inner heliosphere, a large fraction of the velocity profile within MCs is close to a linear function of time. This is 
indicative of a self-similar expansion and a MC size that locally follows a power-law of the solar distance with an exponent called f . 
We derive the value of ( from the in situ velocity data. 

Results. We analyze separately the non-perturbed MCs (cases showing a linear velocity profile almost for the full event), and perturbed 
MCs (cases showing a strongly distorted velocity profile). We find that non-perturbed MCs expand with a similar non-dimensional 
expansion rate (f = 1.05 ± 0.34), i.e. slightly faster than at the solar distance and in the inner heliosphere (£ = 0.91 ± 0.23). The subset 
of perturbed MCs expands, as in the inner heliosphere, at a significantly lower rate and with a larger dispersion (f = 0.28 ± 0.52) 
as expected from the temporal evolution found in numerical simulations. This local measure of the expansion also agrees with the 
distribution with distance of MC size, mean magnetic field, and plasma parameters. The MCs interacting with a strong field region, 
e.g. another MC, have the most variable expansion rate (ranging from compression to over-expansion). 



Key words. Sun: magnetic fields, Magnetohydrodynamics (MHD), Sun: 
Interplanetary medium 



coronal mass ejections (CMEs), Sun: solar wind, 



1. Introduction 

Magnetic clouds (MCs) are a subset of interplanetary coronal 
mass ejections, that show to an observer at rest in the helio- 
sphere: the large and coherent rotation of the magnetic field 
vector, a magnetic field intensity higher than in their surround- 
ings, and a low proton temperature (e.g jKlein & Bur laga 1982; 
Burlaga 1995). They transport huge helical magnetic structures 
into the heliosphere, and are believed to be the most effi cient 
mechanism for re leasing magnetic helicity from the Sun (iRustl 
1994; lLowHl997l) . They are also the most geoeffective events 
to erupt from the Sun, and are associated with observations of 
strong flux decreases for galactic cosmic rays (ITsurutani et al] 



1988; Gosling et al. 1991; Webb et al 
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Despite the increasing interest in improving our knowledge 
of MCs in the past few decades, many properties remain un- 
constrained, as such their dynamical evolution as they travel in 
the solar wind (SW). In particular, duri ng their evolution, MCs 
undergo a significa nt expansion (e.g., iKlein & Burlagal fl982t 
Farrugia et al. 1993). Observations of MCs at different heliodis- 
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tances show th at their radial size increases with the distance to 
the Sun (e.g iKumar & RustUl996r.lBothmer & Schwennll 19981: 
iLeitner et al J20071) . More generally, interplanetary coronal mass 
ejections (ICMEs) are also known to undergo a glo bal expan- 
sion in the radial direc tion with heliodistance (e.g., Wan g et all 
2005b]; iLiu et al.ll2.005h . The expansion properties of MCs and 
ICMEs then differ significantly from the stationary Parker's so- 
lar wind. 



From in situ observations of MCs, the plasma velocity pro- 
file typically decreases as the MC passes through the spacecraft, 
showing a faster speed at the beginning (when the spacecraft is 
going into the cloud, Vr,;,,) and a slower speed at the end (when 
the spacecraft is going out o f the cloud, V« iOUt , e.g., Leppin g et alJ 
2003; Gul isano et alj201Q|) . Moreover, the velocity profile inside 
MCs typically shows a linear profile, which is expected for a 
self-similar expansion (Farrugia et al. 1119931 IShimazu & VandasI 
120021; IDemoulin etail ho08t IDemoulin & Dassol l2009al) . The 
time variation of the size of the full MC can be approximated 
with AVmc = V«,out - Vu; m . However, a parcel of fluid a factor / 
times smaller than the MC would typically have a difference in 
velocity AV = / AVmc (because of the linear profile). 
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Then, as AV is size dependent, it is not the intrinsic expan- 
sion rate of the parcels of fluid. Moreover, the regions near the 
MC boundaries typically undergo the strongest perturbations, so 
that AV is not a reliable measurement of the flux-rope expansion 
rate. 

To measure with greater a ccuracy the expansion of plasma 
in MCs, Dem oulin et al.l (120081) introduced a non-dimensional 
expansion rate (called (, see Sect. 14. Il l, which is computed from 
the local insitu measurements. It was shown that, if ( is constant 
in time, then £ also determines the exponent of the self-similar 
expansion, i.e. the size of a parcel of fluid as r(D) ~ ro/X, with 
D the distance to the Sun. 

From the analysis of the non-dimensional expansion rate (£) 
for a large sample of MCs obs erved in the inner heli osphere by 
the spacecrafts Helios 1 and 2, Gulisa no et al.l (1201 Oh found that 
there are two populations of MCs, one that has a linear profile 
of plasma velocity throughout almost the entire passage of the 
cloud by the spacecraft (the non-perturbed subset) and another 
population that has a strongly distorted velocity profile (the per- 
turbed subset). The presence of the second population was in- 
terpreted as a consequence of the interaction of a fraction of the 
MCs with fast solar wind streams. 

In the present paper we extended the study of Guli sano et al.l 
d2010l) to the outer heliosphere from 1 to 5 AU, analyzing ob- 
servations made by the Ulysses spacecraft. We first describe the 
data used, the method to define the boundaries, the axis direction 
of the flux ropes, and the main characteristics of the velocity pro- 
file (Sect. |2|i. Next, we study statistically the dependence of the 
main MC properties on the solar distance (Sect. |3}- Then, we 
analyze the expansion properties of MCs separately for the per- 
turbed and non-perturbed cases (Sect. 3). The MCs in interaction 
with a strong external magnetic field (such as another MC) are 
analyzed separately case by case in Sect. [5] Finally, we summa- 
rize our results and conclude (Sect. |6). 

2. Data and method 

2.1. Selected MCs 

We derive the generic expansion properties of MCs in the outer 
heliosphere. We do not select MCs based on their physical prop- 
erties (e.g. depending on whether they drive a front shock for 
instance), but rather study a set, as complete as possible, of MCs 
in a given period of t ime. W e select the MCs previously analyzed 
bv iRodriguez et all (2004) based on their thermal and energetic 
properties. This list of 40 MCs covers the time interval from 
February 1992 to August 2002. During the analysis of Ulysses 
data, we find six extra MCs during this time interval, which are 
included at the end of Table I ATI 

Ou r present study has a different scope from Rodri guez et al.l 
(2004) as our main aim is to quantify the expansion properties of 
MCs. We use data from the S olar Wind Observ ations Over the 
Poles of the Sun (SWOOPS, iBame etalJll992l) of plasma ob- 
servations with a temporal cadence of 4 minutes. The magnetic 
field data are from the Vector Helium Magnetometer (VHM, 
Balogh et dll992h . and have a temporal cadence of one second. 

2.2. The MC frame 

The magnetic and velocity fields observations are in the radial 
tangential normal (RTN) system of reference (R, T, N), where 
R points from the Sun to the spacecraft, T is the cross product 
of the Sun's rotation vector with R, and N completes the right- 
handed system (e.g., iFranz & Harper 2002). 



Another very useful system of coordinates is one attached 
to the local direction of the flux rope. In this MC local frame, 
Zcioud is aligned along the cloud axis (with Becloud > in the 
MC central part), and found by applying the minimum variance 
(MV) technique to the normalized time series of the observed 
magnetic field (e.g. Guli sano et al.l200"7l and references therein). 
We define Y c i 0U( j in the direction R x Z c i 0U( j, and X c i ou( j completes 
the right-handed orthonormal base (X c i ou d, Y c i ou d, Z c i ou( j). 

In the MC frame, the magnetic field components have a typ- 
ical behavior that we now describe. The axial field component, 
Xcioud, typically reaches a maximum in the central part of the 
MC, and declines toward the borders to small, or even nega- 
tive, values. The component B y ,cioud typically changes of sign 
in the MC central part. When large fluctuations are present in 
some MCs, then B y , c ioud has an odd number of reversals (e.g. 
ISteed et aDl201 ll) . Finally, the component v ,cioud has its weak- 
est value when the spacecraft crosses the central part of the MC; 
its relativ e magnitude provides an indication of the imp act pa- 
rameter (Gulisa no et al.ll2007t iDemoulin & D asso 2009a). 

Another reason for being interested in the MC frame is that it 
permits us to relate the inward and outward limits of t he crossed 
flux r ope to the conservation of the azimuthal flux (Dass oet al.l 
2006). We neglect the evolution of the magnetic field during 
the spacecraft crossing period (since the elapse d time is small 
comp ared to the transit time from the Sun, see De moulin et al.l 
2008 for a justification), and define the acc umulative flux per 
unit length L along the MC axial direction as dDasso et al . 2007) 

F (?) f ' 

-j- L = J By^uJt') V x , cloud (t') dt' , (1) 

where t m is the time of the front boundary. This time is selected 
as a reference since the front boundary is usually well-defined, 
although any other reference time can be used without any ef- 
fect on the following conclusions. Since B v , c / i«/ undergoes its 
principal change of sign inside a MC, F y (t), then has a main 
extremum and defines the time called f cen tei, which is approxi- 
mately the time of closest approach to the flux rope center. With 
the hypothesis of a local symmetry of translation along the main 
axis, the flux surface passing at the position of the spacecraft at 
time 1 1 is wrapped around the flux rope axis, and is observed at 
time ?2 defined by 

F y {h) = F y (h), (2) 

from the conservation of the azimuthal flux. Within the flux rope, 
this equation relates any inward time (time period for inward 
motion) to its conjugate in the outward (the same but for outward 
motion). 

2.3. Definition of the MC boundaries 

As a first approximation, we choose the MC boundaries from 
the magnetic field in the RTN coordinate system, taking into 
account the magnetic field properties (see Section [TJ. We also 
consider the measured proton temperature and compared them 
to the expected tempera ture in the SW with the same speed 
dLopez & FreemanTl 986) and to S W temperature acc ounting for 
its dependence on solar distance (Wan g et aLll200 5b). However, 
a lower temperature than expected is only used as a guide, since 
we found that the magnetic field has sharper changes than the 
temperature (observed in the time series of data), and so the mag- 
netic field defines more precisely the flux rope boundaries. 

The front, or inward, boundary of a MC is typically the easi- 
est to determine. It is the limit between the fluctuating magnetic 
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Fig. 1. Examples of two MCs with velocity profiles that are not significantly perturbed (i.e. that have an almost linear dependence on 
time), (a) MC 13 travels at low latitude (2°) at D = 5.3 AU in a slow SW and (b) MC 30 travels at high latitude (-66°) at D = 3 AU 
in a relatively slow SW for such high latitude. The vertical dashed lines define the MC (flux rope) boundaries, and the vertical dotted 
line defines the rear boundary of the MC back (Sect. 12.31 a back is only present on the right panels). The three top panels show the 
magnetic field norm and its two main components in the MC frame computed with the MV method (Sect. 12.21 ). Z? z ,mv is the MC 
axial magnetic field component. B y ,mv is the magnetic field component both orthogonal to the MC axis and to the radial direction 
from the Sun (R). The solid red line represents F y , which is the accumulated flux of B y MV [Eq. ©]. The bottom panel shows the 
observed velocity component in the radial direction (R). A linear least squares fit of the velocity (green line) is applied in the time 
interval where an almost linear trend is present within the MC. 




Fig. 2. Examples of two MCs with perturbed velocity (i.e. deviating from a linear function of time), (a) MC 1 1 travels at low latitude 
(6°) at D = 5.2 AU in a slow SW and (b) MC 36 travels at high latitude (75°) at D = 2.3 AU in a fast SW. We use the same quantities 
and drawing convention as Fig.Q] 



field of the sheath and the smoother field variation within the 
flux rope. A current sheet is typically present between these two 
regions, and the front boundary is defined to be there, at a time 
denoted t m . 

The rear, or outward, boundary is more difficult to define, 
particularly because many MCs have a back region with phys- 
ical properties that are in-between those of the MC and SW 
hence the transition between the flux rope and the perturbed 
SW i s frequentl y not as well defined as at the inward bound- 
ary. iDasso et al .] (120061 [2007 ) concluded that this back region 



is formed, during the transit from the Sun, by reconnection of 
the flux rope with the overtaken magnetic field. We use the con- 
servation of the azimuthal flux between the inward and outward 
boundary, Eq. ©, to define the outward boundary. When a back 
region is present, the time f out of the outward boundary is de- 
fined as the first solution of F Y (t out ) = F y (t m ) = 0. However, this 
condition depends on the determination of the MC frame thus an 
iterative procedure is needed, as follows. 

We start with approximate MC boundaries defined in the 
RTN system, then perform a MV analysis to find the local frame 
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of the MC. Next, we analyze the magnetic field components in 
the local frame, and redefine the boundaries according to the az- 
imuthal flux conservation expected in a flux rope [Eq. ©J. In 
most cases, the temporal variations in the magnetic components 
indicate that the rear boundary needs to be change its location. 
The MV analysis is then performed only on the flux rope inter- 
val. This is an important step because an excess of magnetic flux 
on one side of the flux rope typically introduces a systematic 
bias in the MV directions. Finally, we perform the same proce- 
dure iteratively until convergence is achieved. This defines more 
reliably determined boundaries and orientation of the flux rope. 

As at the front boundary, we expect to find a current sheet at 
the outward boundary (since it separates two regions with differ- 
ent origins and physical properties). A current sheet is typically 
found at about time f out in most analyzed MCs . The bending of 
the flux rope axis is one plausible cause of the small difference 
in time between f out and the closest current sheet. 

In the outward branch of a few MCs, B VjC ioud has an impor- 
tant current sheet and later on it fluctuates before F y (f) finally 
vanishes. This indicates that there is a deficit of azimuthal flux 
in the outward branch compared to the inward one. In this case, 
t out is defined by the observed time of the current sheet, and the 
flux rope extension is in the time interval [fi n ,f utL where t m is 
redefined by F y (t m ) = F y {t out ). The iterative procedure described 
above is then applied to improve f; n and the flux rope orientation. 

2.4. MC groups 

Following Gulisan o et al.l d2010h . we divide the data set into 
three groups: non-perturbed MCs, perturbed MCs, and MCs in 
interaction (noted interacting). 

The interacting group contains the MCs that have a strong 
magnetic field nearby (e.g. another MC). From MHD sim- 
ulations, their physical properties are exp ected to differ sig - 
nificantly from non-interacting MCs (e.g.. IXiong et al.l 12007). 
Observations, indeed confirm that MCs in the process of inter- 
actin g could have a different behavior from non-interacting MCs 
(e.g. lWang et d1l2005atlDasso et al.ll2009h . 

The MCs in the non-perturbed and perturbed groups have 
less disturbed surroundings than MCs in the interacting group. 
The MCs in the non-perturbed group have an almost linear ve- 
locity profile in more than 75% of the flux rope, while in the 
perturbed MCs this condition is not satisfied. Two examples of 
each group are shown in Figs. 1 1121 the label (number), the group, 
and the main properties of each MC are given in Table IA.fl 

3. Magnetic cloud size and mean properties 

We now analyze the size and mean properties, along the space- 
craft trajectory, of MCs as a function of solar distance D. 
Because of the trajectory of Ulysses, a correlation between D 
and the absolute value of the latitude is present in the ana- 
lyzed data (Fig. |3j. However, the results of the analysis that 
we now describe do not change significantly when MCs are 
grouped in terms of the latitude, so that this correlation between 
D and latitude does not affect our conclusions. All the studied 
properties are affected by the expan sion achieved from the Sun 
to Ulysses. From models (see e. g., IChenl fl996: Ku mar & Rust! 
1996; Demoulin & Dasso 2009a), all properties are expected to 
have a power-law dependence on D, we then perform a linear fit 
for each property plotted on a log-log scale (Figs. 14151 . for the 
non-perturbed and perturbed groups. Since MCs in interaction 
are strongly case-dependent, their properties are not compared 
(fitted) with the global models here. 
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Fig. 3. Spatial localization of the MCs listed in Table I A. II The 
MCs are separated into three groups: non-perturbed (filled green 
circles), perturbed (empty blue squares), and in the process of 
interaction (red stars). 



3.1. Magnetic cloud size 

With the MC boundaries defined in Sect. I2.3I we compute the 
size S in the direction R as the product of the duration of the MC 
and its velocity at closest approach to the center of the flux rope. 
Compared to most previous works (see Sect. [1}, this typically 
defines a smaller size as the back region is not included. 

Since the majority of previous studies include all MCs in- 
dependently of the groups we defined above, we first compare 
them to our results with all MCs included (Table[T]i. In the range 
[1.4,5.4] AU, we find an exponent that is lower (0.49+0.26) than 
previous studies on MCs or ICMEs in the inner and outer helio- 
sphere (lower part of Table [T). Taking into account only non- 
perturbed and perturbed MCs, the exponent is slightly higher 
(0.56 + 0.34), but still lower than previous studies, including our 
previous result obtained with Helios data (0.78 + 0.12). 

Next, we analyze separately the non-perturbed and perturbed 
MCs. In the range [1.4,5.4] AU, we find comparable exponents 
to our previous results obtained with Helios data for both non- 
perturbed and perturbed MCs (Table [T). Perturbed MCs have a 
significantly smaller size on average than non-perturbed MCs 
and this tendency increases with solar distance (Fig.|4^). We find 
that perturbed MCs typically have smaller sizes S , by a factor 
a 1.3, than non-perturbed MCs. 

The MCs in interaction (interacting MCs) have typically 
smaller sizes S than the non-perturbed MCs (Fig. |4^). We in- 
terpret this result as a direct effect of the compression and/or 
reconnection with the interacting MC or strong magnetic-field 
structure. 

3.2. Magnetic field strength 

We define the average field < B > within the flux rope bound- 
aries defined in Sect. l2.3l and proceed as above with S . The value 
of < B > decreases far more slowly with distance in the Ulysses 
range (-1.18 + 0.27) than the Helios range (-1.85 + 0.07) for 
the non-perturbed and perturbed MCs (a similar result is ob- 
tained when also including the interacting MCs). Finally, com- 
bining our Helios and Ulysses results, we have an intermediate 
exponent (-1-5 + .1). T his last re s ult is closer to the results of 
Wan g et al.l (|2005b) and ILiu et al.1 (12005b obtained for a larger 
set of ICMEs across the same range of distances (bottom of 
Tabled}. 
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Table 1. Exponents e in the fit with a power law D e for S, < B >, 
and < N p > (found by performing a linear fit in log-log plots, 
Fig.©. 



range of D a 


S b 


< B > b 


<N p > b 


Ref. c 


all MCs 










[1.4,5.4] 


0.49±0.26 


-1.20+0.20 


-1.70+0.34 


□ 


non-perturbed and 










perturbed MCs 








□ 


[0.3, 1] 


0.78±0.12 


-1.85+0.07 




[1.4,5.4] 


0.56+0.34 


-1.18+0.27 


-1.70+0.43 


□ 


non-perturbed MCs 










[0.3, 1] 


0.89±0.15 


-1.85+0.11 






[1.4,5.4] 


0.79+0.46 


-1.39+0.92 


-2.24+0.66 


□ 


perturbed MCs 










[0.3, 1] 


0.45+0.16 


-1.89+0.10 




111 


[1.4,5.4] 


0.54+0.48 


-1.14+0.26 


-1.40+0.50 


□ 


MCs 










[0.3, 1] 






-2.40+0.30 


® 


[0.3, 1] 


1.14+0.44 


-1.64+0.40 


-2.44+0.46 




[0.3,4.] 


0.97+0.10 


-1.80 


-2.8 


® 


[0.3,4.8] 


0.78+0.10 






® 


[0.3,5.4] 


0.61+0.09 


-1.30+0.09 


-2.62+0.07 


® 


[1.4,5.4] 




-0.88+0.22 




a 


ICMEs 










[0.3,5.4] 


0.61 


-1.52 




m 


[0.3,5.4] 


0.92+0.07 


-1.40+0.08 


-2.32+0.07 


m 



[3l iBothmer & Schwenn] dl998h. HI iLeitner et al.l (120071). 
iKumar & Rusj d!996h . ISl lwang et alj T2005b). ITl lLiu et alN2005h . 



The non-perturbed MCs have a slightly steeper decrease in 
magnetic field strength with distance than the perturbed MCs 
(Fig. H]}). This result is coherent with the conservation of mag- 
netic flux combined with a slightly larger increase in size with 
distance of the non-perturbed MCs compared to the perturbed 
ones. We also find that the perturbed MCs typically have a 
stronger magnetic field, by a factor of between « 1 .4 and 1 .7, 
than non-perturbed MCs (Fig. |4j>), a result that is consistent 
with the lower expansion rate for perturbed MCs (see Sect. 13. II 
andl43T>. 

The interaction has a weaker effect on < B > than on the 
size S : There is only a weak tendency for interacting MCs to 
have stronger < B > than non-perturbed MCs (Fig.|4j)). 
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log D 
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' d The range of solar distances is in astronomical units (AU). In the 
middle part of the table, MCs are separated in non-perturbed and 
perturbed MCs. 

b The exponents, e, are given for the size 5 in the R direction, the 
mean magnetic field strength < B >, and the mean proton density < 
N p >. All quantities are computed within the flux rope boundaries. 
c We compare our results with our previous results in the in- 
ner heliosphere and with other studies of MCs and ICMEs. Th e 
results are from: [TJ pre s ent w ork, l2l iGulisa no et al.j (| 2010t) . 



2.5 



1 .5 



0.5 



-0.5 



-1 .5 



log <B>=~1 .39 log D+0.94 
log <B>=-1 .1 4 log D+1 .01 
log <B>=-1.18 log D+0.94 



r=0.66 
r=0.76 
r=0.65 



(non perturbed) 

(perturbed) 

(both) 



(b) 




0.4 
log D 



0.7 0.8 



log <N 
log <N p 
log <N 



-2.24 log D+0.63; 
-1 .42 log D+0.43; 
-1 .70 log D+0.49; 



=0.70 
=0.57 
=0.60 



(non perturbed) 

(perturbed) 

(both) 



(C) 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

log D 

Fig. 4. Dependence on the solar distance, D in AU, of (a) the 
radial size in the R direction (AU), (b) the mean magnetic field 
strength (nT), and (c) the mean proton density (cm -3 ) with log- 
log plots. The MCs are separated in three groups: non-perturbed 
(filled green circles), perturbed (empty blue squares), and in in- 
teraction (red stars). The straight lines are the result of a least 
squares fit for non-perturbed (thin continuous green line), per- 
turbed (dashed blue line), and for non-perturbed and perturbed 
MCs (thick continuous black line). The fits and the absolute 
value of the correlation coefficients, r, are added at the top of 
each panel. 



3.3. Proton density 

We define the average proton density < jV p > as above for the 
magnetic field strength, i.e. within the flux rope boundaries. For 
the non-perturbed and perturbed MCs, the density decreases sig- 
nificantly less rapidly with distance (exponent = -1.70 + 0.43) 
than in previous studies (exponents in the range [-2.8, -2.3]) for 
both the inner heliosphere and a combination of both the inner 
and outer heliosphere. However, we are unaware of results for 



the outer heliosphere alone (apart Fig. 10 of iLeitner et al.l l2007. 
but no fit is provided). 

Nevertheless, for the non-perturbed MCs the density de- 
creases with D following a powerlaw that is in closer agreement 
with previous studies. In contrast, for perturbed MCs the density 
decreases much more slowly with distance. This agrees quali- 
tatively with a lower expansion of S and a lower decrease in 
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Table 2. Average and dispersion in the main parameters of MCs from our present study and lDu etal 1 (120101) . 



Group a 


1/ b 


S IAV C 


< B >l AU c 


< N p >iau c 


<P P > 


(i 

7 


i c 




km/s 


AU 


nT 


cm -3 








all MCs 


470±110 


0.16±0.12 


12+7 


9+8 


0.19+0.17 


68+17 


0.77+0.75 


non-perturbed and perturbed 


480±120 


0.18±0.14 


14+8 


9+8 


0.18+0.19 


68+20 


0.63+0.59 


non-perturbed 


480±130 


0.24±0.16 


11+9 


7+10 


0.17+0.18 


63+24 


1.05+0.34 


perturbed 


480±110 


0.13+0.09 


16+6 


11+7 


0.19+0.20 


72+15 


0.28+0.52 


interacting 


470+130 


0.12+0.05 


10+6 


9+9 


0.22+0.14 


66+13 


1.00+0.93 


MC f 


500+100 


0.26+0.19 


11+6 








0.59+0.51 


non-MC ICME f 


480+ 90 


0.22+0.21 


5+4 








0.68+0.48 



il The MCs are separated into three groups: non-perturbed, perturbed, and interacting. 
b V c is the velocity at the closest distance from the MC center. 

G The size and the mean field are normalized to 1 AU according to Eq. ([8) and the proton density with < N p >i au =< N p > D 12 , where D is the 

solar distance in AU. 
d y is the angle between the MC axis and the radial direction (R). 
e f is the nondimensional expa nsion rate [Eq. d5}]. 
f Computed from the results of iDu et all ( 120101) 7 
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Fig. 5. Dependence on the solar distance, D in AU, of (a) the 
mean proton temperature < T p > (K) and (b) the mean proton B. 
The drawing convention is the same than in Fig. [4] 



< B > with D, while the difference between non-perturbed and 
perturbed MCs is more significant for the average proton density. 
Even if there is qualitative agreement, the systematically smaller 
decrease in < N p > with D than in previous studies could be re- 
lated to the substant ial variability of iVp in MCs (e.g. at 1 AU see 
iLepping et al.ll2003l) so that < N p > is case dependent and could 
be biased with our relatively small number of MCs. When the 
same MC is observed at two distances, this bias is removed. For 
MC 15, we computed from the observed data a steeper slope of 



2.9 by using the mean o bserved density at both ACE and Ulysses 
(Nak wackietai1l201lh . 

The effect of the interaction on < N p > is weak, and com- 
parable to the effect on < B >: there is only a weak tendency 
for interacting MCs to have stronger < N p > than non-perturbed 
MCs (Fig. [4};). While there are too few interacting MCs at low 
D to enable us to draw a strong conclusion, the trend of a global 
decrease in < A^ p > with D predominates over the effect (typi- 
cally compression) of an interaction with another MC or a strong 
B region. 



3.4. Proton temperature 

The average proton temperature decreases only weakly with dis- 
tance (Fig. The temperatures of the non-perturbed and per- 
turbed MCs have dependences close to D~~° M and D~ 0M , respec- 
tively. This contrasts wit h the decrease found in the SW, which 
is typ ically around D~ 0J (iGazis et al.ll2006h iRichardson & Canel 
119951) . Moreover, the interacting MCs are not significantly hot- 
ter than non-perturbed and perturbed MCs (apart from two hot 
interacting MCs, Fig.|5^). 

The above result also dilfers from the dec r ease in temper- 
ature given by « D~ 1 - 6 found by iLeitner et al.l (120071) for MCs 
observed in th e range r0-3, 5.41 A U. However, our results agree 
with those of Leitn er et al.l (|2007) since their Fig. 1 1 shows no 
decrease in proton temperature with distance for MCs observed 
by Ulysses. We conclude, in agreement with previous results, 
that the proton temperature of MCs in the outer heliosphere de- 
creases with D much more slowly than in the inner heliosphere 
and than in the SW. 

Moreover, MCs have a much larger volume expansion than 
the SW, the main difference being that MCs have a significant 
expansion in the radial direction (Table [TJ and Fig. [6] below) 
in contrast to the SW, which has no mean radial expansion. 
Considering an adiabatic expansion implies that MCs should be- 
come even cooler than the SW with distance. The opposite is in 
fact observed, which implies that there is a far greater amount of 
heating in the MCs than the S W. A significant source of this heat- 
ing could be a turbulent cascade to small scales. Since B p « 1 
in MCs (see Sect. 13.51 ). the dissipation of a small fraction of the 
magnetic energy provides a large increase in the plasma temper- 
ature, while the same amount of magnetic energy in the SW pro- 
vides only a small increase in the plasma temperature (B p « 1). 
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3.5. Proton plasma B 

The proton plasma B, denoted /? p , is defined as the ratio of 
the proton thermal pressure to the magnetic pressure. Thus, if 
the decay of the proton density, proton temperature, and mag- 
netic field intensity are true power laws (respectively, N p oc 
£>"», T p oc D " T , and B oc D _ " B ), the decay of B v will be 
B p oc D- n "-" T+2 " B , which in our Ulysses study corresponds to 
an expected increasing function with the heliodistance, such as 



oc D 



■1.7-0.31+2*1.18 



Z)+u.j5 for the combined sample of 



non-perturbed and perturbed MCs (Figs. Ht,c, |5^). Another es- 
timation can be realized by supposing an isotropic expansion of 
MCs with distances that increase as D m . Assuming the conser- 
vations of both mass and magnetic flux (i.e., an ideal regime), 
we theoretically expect an evolution such that B p oc £)-3m-n r +4m j 
which is Bp oc D +0 25 for non-perturbed and perturbed MCs and 
B p oc D +0 - 45 for non-perturbed MCs, and then comparable to the 
previous estimation. 

We compute the mean value of B p within each studied cloud, 
finding that < B p > is lower than unity for all MCs in our sam- 
ple, and has a tendency to slightly increase with solar distance 
for non-perturbed and perturbed MCs. We find that B p oc D° 35 
for the combined groups of MCs (Fig. \5&), so a very similar de- 
pe ndence to that f ound above assuming exact power laws. 

iLeitner et al.l (120071) found a slight decrease, as D~ 0A , for 
MCs in the range [0.3,5.4] AU. However, this global tendency 
is not present for MCs observed with Ulysses, as their Fig. 12 
shows a tendency of increasing /3 P with D. We then conclude 
that Bp has a different dependence on D in both the inner and 
outer heliosphere. 

Next, we analyze separately the dependence of B p on D for 
each MC group. The value of /3 p is almost independent of D for 
the non-perturbed MCs, and the above increasing trend with dis- 
tance can be mainly attributed to the perturbed MCs. The MCs 
that are interacting typically have a larger B p at all D (Fig.|5j3). 

4. Expansion rate of MCs 

Magnetic clouds typically have a velocity profile close to a lin- 
ear function of time with a larger velocity in the front than at 
the rear. Hence MCs are expanding magnetic structures as they 
move away from the Sun. In this section, we characterize their 
expansion rate. 

4.1. Non-dimensional expansion rate 

The measured temporal profile V R (t) of each MC is fitted using 
a least squares fit with a linear function of time 



V S M(f) = V oM + (.dV R /dt) m t, 



(3) 



where (dVs/df)g t is the fitted slope. We always keep the fitting 
range within the flux rope, and restricts it to the most linear part 
of the observed profile. The observed linear profile of Vr indi- 
cates that different parts of the flux rope expand at about the 
same rate in the direction R. 

The linear fit is used to define the velocities V/j,fi t (f/„) and 
VRfi t (t „r) at the flux rope boundaries (Sect. |23T >. We then define 
the full expansion velocity of a flux rope as 



AV R = V RM (t in ) - V RM (t out ) . 



(4) 



For non-perturbed MCs, AVr, is very close to the observed ve- 
locity difference V«(f,„) - Vr(W)> see e.g. the lower panels of 
Fig-Q] For perturbed MCs, this procedure minimizes the effects 



of the perturbations within the flux rope and especially close to 

the MC bound aries. 

Fo llowing iDemoulin & Dassol (l2009al) and iGulisano et al.l 
(2010), we define the non-dimensional expansion rate as 



AVr D 

~V} 



(5) 



which defines a scaling law for the size S of the flux rope (along 
R) with the distance (D) from the Sun as S oc D ( . This simple 
interpretation of ( is obtain ed with the following two si mplifica- 
tions which were justified in Dem oulin & DassoT(l2009al) . Firstly, 
we neglect the aging effect (the front is observed at an earlier 
time than the rear, so when observed at the front the flux rope is 
smaller than when it is observed at the rear), and secondly, ( is 
approximately constant during the time interval of the observed 
flux rope. 

While £ is a local measurement within the flux rope, it pro- 
vides a measure of the exponent m of the flux rope size if it 
follows a power law 



taking the temporal derivative of S , we find that 

AtV}. 



dS dS dD S 

— = AV R = « m—V c 

dt dD dt D 



D 



(6) 



(7) 



We then find a relation that is equivalent to Eq. (0, so for self- 
similarly expanding MCs we have m « £. 

4.2. Expansion of non-perturbed MCs 

The main driver of M Cs expansion was identified by 
Dem oulin & Dassol (l2009al) as the rapid decrease in the total S W 
pressure with solar distance . They followed the force-free evolu- 
tion away from the Sun of flux ropes with a variety of magnetic 
field profiles and assuming either ideal MHD or fully resistive 
relaxation that preserves magnetic helicity. Within this theoreti- 
cal framework, they demonstrated that a force-free flux rope has 
an almost self-similar expansion, so a velocity profile is almost 
linear with time as observed by a spacecraft crossing a MC (e.g. 
Figs. 1112b . In the case of a total SW pressure behaving as D~ np , 
they also found that the normalized expansion rate is t, « np/4. 
These results apply to a progressive evolution of a flux rope in 
a quiet SW, so the non-perturbed MCs are expected to have the 
closest properties to these theoretical results. 

The mean value of £ for non-perturbed MCs, observed in 
the range [1.4,5.4] AU, is < £ >« 1.05 + 0.34. This is slightly 
above the mean values, * 0.80 ± 0.18 and « 0.91 + 0.23, 
found at 1 AU and in the range [0.3,1 .1 AU, respectively 
( Demoulin et al .1 120081: IGulisano et alJ[2010h . This is a small in- 
crease compared to the change in D, which is a factor in the 
defini tion of ( [Eq. © I. Since V c has no significant dependence 
on D (IWang et al.l l2005b). the increase in D is compensated for 
mainly by the decrease in AVr/ At. 

The slightly higher mean < £ > in the range [1.4,5.4] AU 
is mostly due to the higher £ values found at smaller D values 
(Fig. [6^), so in the region that most closely matches that of both 
previous studies. Owing to the characteristics of the Ulysses or- 
bit, there are indeed only a few observed MCs for D < 3 AU 
(Fig. ID so that the decrease with D in the linear fit of f could be 
due to the specific properties of the few detected MCs at these 
lower distances. We also recall that D is strongly correlated with 
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Fig. 6. The panels (a-d) show the correlation analysis that tests for the dependence of the non-dimensional expansion factor £ 
[Eq. (0] on other MC parameters. The drawing convention is the same as in Fig. [4] Both S i au and < B >i au are normalized to 1 
AU using the size and field strength dependence on the distance, according to the relationship given in Figure [4] for non-perturbed 
MCs, see Eq. (HI). 



the heliolatitude \6\, Fig. [3] so that the dependence ((D) shown in 
Fig. [(J could also be an effect of latitude. Nevertheless, we em- 
phasize that the < £ > found here with Ulysses is only slightly 
higher than the values found previously at smaller distances from 
the Sun. These differences are small compared to the signifi- 
cantly different SW properties of the slow and fast SWs pre- 
dominantly present at low and high latitude, respectively. 

As in the previous study, we test the possible dependence 
of ( on the MC properties. Both the MC size and field strength 
depend strongly on D (Fig. |4^,b). We use the fits found for the 
non-perturbed MCs to remove, on average, the evolution with D. 
Hence we define values at a giving solar distance, here taken at 
1 AU, with the relations 



Si 
< B >, 



-0.79 



- S D 

- < B > D 



-1.39 



(8) 

The mean MC spee d depends only we akly on D in the outer he- 
liosphere (see e.g. iWang et al. 2005b, for a large set of ICMEs, 
including MCs), hence we use below the measured velocity of 
the center V c . 

For non-perturbed MCs we find that ( is almost indepen- 
dent of < B >iau (Fig, [(pi), as found in the r ange [0.3, 1.] AU 
(compared to Fig. 4d of iGulisano et al.ll2010l) . A small differ- 
ence is that f decreases slightly with V c and S \ au (Fig. Et>,c), 
while a nearly constant value was found in the ran ge [0.3, 1.] 
AU (compared to Fig. 4b,c of IGulisano et al. 2010). This dif- 
ference is probably not due to a latitude dependence since the 
MCs observed at |f?| < 25° or > 25° are evenly distributed in the 
three plots when we analyze the two groups independently [not 



shown], hence similar results are obtained. Finally, using alter- 
native exponents in Eq. (O, in the range given in TableQ] induces 
only slight changes in the linear fits shown in Fig. [6] particularly 
in the magnitude of the slope. 

4.3. Expansion of perturbed MCs 

The main difference between perturbed and non-perturbed MCs, 
is that the former has a lower expansion rate and a larger dis- 
persion, < ( >= 0.28 ± 0.52, compared to < £ >= 1.05 ± 0.34 
for non-perturbed MCs (see also Fig. |7];,d). These results are 
com parable to the results found in the distance range [0.3, 1] AU 
bv IGulisano et al.1 d2010h since they found < £ >= 0.48 + 0.79 
and < £ >= 0.91 + 0.23 for perturbed and non-perturbed 
MCs, respectively (Fig.|7^,b). This lower expansion rate for per- 
turbed MCs is indeed expected from MHD simulations if the 
origin of the perturbation is an overtaking magnetized plasma, 
which compresses t he MC, hence decreases its expansion rate 
(IXiong et al.ll2006allbh . 

A substantial overtaking flow, with a velocity larger than at 
the front of the MC, is found in 12 of the 20 perturbed MCs. 
These strong flows are located behind the observed MC, and 
typically a shock, or at least a rapid variation in V, is present 
at the rear boundary of the back region, as illustrated by MC 1 1 
in Fig.|2^. This is the main difference from MCs observed in the 
range [0.3, 1] AU, as the overtaking flow typically entered deeply 
into the perturbed MCs, and moreover, a strong shock was typ- 
ically present within the pe rturbed MC (compare Fig. |3£ with 
Fig. 2 IGulisano et a 1. 2010). This difference can be interpreted 



8 



Gulisano et at : Expansion of magnetic clouds in the outer heliosphere 



Helios 1 & 2 
10- Non-Perturbed MCs 




Helios 1 & 2 
io- Perturbed MCs 




Ulysses 

Non-perturbed MCs 




Ulysses 

Interacting MCs 




(eh 











60 


- MC analyzed 


(f) 


40 


by Du et al. 






20 
















60 - ICMEs analyzed 
40 _ by Du et al. 

20 - 



(g) 



0^ 
-2 



-1 



Fig. 7. Histograms comparing £ of MCs and non-MC IC MEs. 
(a,b) MCs from Helios spacecraft dGulisano et al.ll20Toh . (c-e) 
MCs from Ulysses spacecraft (present work). (f,g) MCs and non- 
MC ICMEs from Ulysses with app roximate ( values computed 
from the results of lDu et al.l(l2010l) . 



with the MHD simulations of lXiong et aD (l2006allbh . as follows. 
The observations at Helios distances correspond to the interact- 
ing stage wh en the overtaking sh ock travels inside the MC (e.g. 
see Fig. 5 of Xion g et al1l2006ah . At larger solar distances, the 
overtaking shock exits from the front of the MC, hence it is not 
observed inside the MC with Ulysses data. However, some of 
the overtaking flow is present behind the MC. The compression 
of the overtaking flow thus persists, and most perturbed MCs are 
under-expanding . 

Additional evidence of a different evolution stage at Ulysses 
compared with Helios is that £ shows no significant dependence 
on D (Fig. |6^), while the mean tendency of f was to increase 
with D and to reach the valu e of < ( > at 1 AU o f unperturbed 
MCs for Helios MCs (Fig. 4a lGulisano etai]l2010h . 

Only one perturbed MC expands slightly faster than the 
mean expansion rate of unperturbed MCs (Fig. |6), while a 
stronger over-expansion is present in some MCs in interaction. 
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Fig. 8. Perturbed and non-perturbed MCs have a remarkably dif- 
ferent behavior of £ when they are plotted as a function of AVr. 
The drawing convention is the same as that in Fig. [4] for panel 
(a), while the perturbed MCs are grouped in terms of velocity 
and distance in panel (b), as indicated by the legend. 



These, over-expanding flux ropes can be identified in the above 
MHD simulations in the late stage of evolution, as follows. The 
overtaking magnetized plasma progressively flows on the MC 
sides and overtakes the MC. When only a weak overtaking flow 
remains at the MC rear, the expansion rate could increase. The 
MC internal magnetic pressure is indeed higher (owing to the 
previously applied compression) than the pressure value for an- 
other MC at the same position and with no overtaking flow be- 
fore. This over-pressure drives a more rapid MC expa nsion (e.g. 
see the summary of two simulations in Fig. 7 of IXiong et al.l 
2006a). Hence, the expansion rate of perturbed MCs depends on 
the interaction stage in which they are observed, (see the cartoon 
in Figure 6 of|Gulisano et al. 2 010h . 

Some properties of £ for perturbed MCs are still similar in 
the inner and outer heliosphere. The value of £ shows an increase 
with both y c and S i au (Fig.|6j3,c), and a decrease with < B >\ au 
(Fig. j 6ji) as found with Helios data (Fig. 4b.c.d lGulisano et"al1 
2010). However, the slope values are smaller by a factor w 2.2, 
and 6, respectively. Finally, we note that the above differences 
are not due to a difference in the range of parameters since V c , 
S i au, and < B >i au have similar ranges for Helios and Ulysses 
MCs. 

As for Helios MCs, £ does not depend on the expansion ve- 
locity AVr [defined by Eq. ©] for non-perturbed MCs observed 
by Ulysses, while f is strongly correl ated with AVr for per turbed 
MCs (compare Fig. [8^ with Fig. 5 in lGulisano et al.ll2010l) . This 
result extends previous Helios results to large distances. Non- 
perturbed MCs indeed have an intrinsic expansion rate £ (which 
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is given by the decrease in the total SW pressure with solar dis- 
tance). This is not the case for perturbed MCs, which are in a 
transient stage, so AVr cannot be com puted by dS /d t as do ne in 
Eq. 0. Following the derivation of Gulisano et al. (2010|), we 
instead find that 



^perturbed 



AV R D AV R D<»D l 



S V r 



S<)V r 



(9) 



where the MC size So is taken at distance Dq. For perturbed 
MCs, one then expects that f is linearly dependent on AVr, 
which is indeed true (Fig. [SK). Interacting MCs are strongly per- 
turbed, and indeed f depends more strongly on AVr than the 
group we call perturbed. 

Moreover, from Eq. (0, a dependence of £ on both D and V c 
is expected. Within the observed ranges, both parameters have a 
comparable effect on the slope. Since we have a too small num- 
ber of MCs to perform a test to assess the dependence on both D 
and V c , we only considered the two groups where D and V c ap- 
pear to have a combined effect on the slopes. Hence, these two 
groups have larger slope differences, and luckily, are also the 
most numerous groups. The first group is defined by V c < 500 
km/s and D > 4 AU, or more precisely < V c >~ 410 km/s and 
< D >x 5.1 AU. The second group is defined by V c > 500 
km/s and D < 4 AU, or more precisely < V c >~ 640 km/s 
and < D >w 2.6 AU. We use a typical So = 0.2 AU value at 
Do = 1 AU, and m = 0.54 as found for perturbed MCs (Table0. 
The first group has an expected slope of 0.026, which is com- 
parable to the one derived from the data (0.028, Fig. 0)). The 
second group is expected to have a smaller slope value of 0.012 
and the observations do indeed give a slope that is even lower 
by a factor 1.5. One can also compare the ratio of predicted to 
observed slopes to eliminate the influence of Sq: the observed 
ratio of slopes is a factor 1 .6 larger than the predicted one, and 
decreases to 1.4 if we use m = 0.28, the mean £ value found for 
perturbed MCs (Table0. Taken into account the uncertainties in 
the MC parameters, these results indeed show that Eq. (0 illus- 
trates that there is a quantitative dependence of £ on both D and 
V c for perturbed MCs. 

4.4. Global radial expansion 

The determination of £ provides the expansion rate of the MC 
at the time of the observations, while its size S depends on the 
past history of (, If S were observed at another solar distance, 
this would provide global information on the expansion between 
the two distances (as well as another piece of local informa- 
tion about Q. However, the observation of the same MC by two 
spacecraft is a rare event since it requires a close alignment of 
the spacecraft positions wi th the propagation direc tion of the MC 
(one case was analyzed bv lNakwacki et al]l201ll) . Nevertheless, 
the statistical evolution of S with D provides a constraint on 
the mean evolution of MCs. The flux rope size of non-perturbed 
MCs have a mean dependence close to (Fig. [4^) 



S ocD' 



,0.79±0.46 



(10) 



If a MC evolves while £ remains constant, this implies that 
S oc ZX (Sect. 14. 11 1. The non-perturbed MCs are expected to have 
the most stable £ value. Their mean £ value, 1 .05 +0.34 (Table0, 
indicates that they are undergoing a more rapid expansion than 
indicated by the global results of Eq. (TTQb . The possible differ- 
ence is that the flux rope size decreases with distance owing to 
a reconnection with the encountered SW magnetic field. The ev- 
idence of this reconnection is a back region observed in a large 



number of MCs (e.g. a back region is present in the three MCs 
shown in Figs. II 121) . If the flux rope size were to decrease by a 
factor a: 1.5 between 1.4 and 5.4 following a power law of D, 
it would decrease S by D° 3 , so sufficiently to explain the above 
difference between the mean evolution of S [Eq. (TTOb l and < £ > 
for non-perturbed MCs. 

An alternative possibility is that, at least, some of the non- 
perturbed MCs were in fact previously perturbed so that their 
f value was lower. This would be indicative of a lower global 
expansion of the MC, hence a weaker dependence of S on D. 
We indeed find that < £ >= 0.77 + 0.75 for non-perturbed, 
perturbed, and interacting MCs taken together, which is much 
closer to the exponent in Eq. ( ITOb than < f > for the unperturbed 
MCs alone. According to this view, the unperturbed, interacting, 
and perturbed classifications are local ones, which are only valid 
around the time of the observation, and a MC changes its group 
(non-perturbed/perturbed/interacting) as it travels away from the 
Sun. Finally, both of the aforementioned possibilities could be 
happening. 



4.5. Magnetic cloud and ICME expansion 

We now analyze a broader set of MCs and ICMEs from the re- 
sult s published bvlDu et all d2010l) . They defined ICMEs follow- 
ing I Wang et all d2005bl) as the time intervals where the proton 
temperature, T p , was less than half the expected temperature in 
the SW with the same speed when measured in Kelvin degrees. 
They found that about 43% of the identified ICMEs are MCs. 
Their Ulysses data set covers the time interval from January 
1991 to February 2008. The mean values of the principal pa- 
rameters of this data set are given in Table [2] Both groups have 
similar characteristics to the MCs studied in this paper, with the 
exception of a significantly lower magnetic-field strength in non- 
MC ICMEs. 

iDu et al.l d2010h performed a linear fit to the radial velocity 
through the full ICME (defined by T p , while we restrict our- 
selves to the interior of the flux rope in our study). They de- 
fined AV as the speed difference between the leading and trail- 
ing edges of the ICME. The d istribution of A V is quite similar 
for MCs and non-MC ICMEs dDu et alj|20lol Fig. 4). The main 
differences are merely that few MCs have a large AV and that 
more non-MC ICMEs have very small AV. 

In their table of ICMEs, IDu et alJ <l2010h provided all the 
quantities to compute ( [Eq. {5)]. We find Gaussian-shaped his- 
tograms of £ for both MCs and non-MC ICMEs (Fig.0,g). They 
are very similar; in particular, their means and dispersions are 
comparable = 0.59 + 0.51 for MCs and £ = 0.68 + 0.48 for 
non-MC ICMEs). This is compatible with a common idea that 
non-MC ICMEs have the same properties as MCs, but are sim- 
ply observed in situ near their periphery in such a way that the 
internal flux rope is missed. 

These results are also comparable to our results when we 
group all MCs together, particularly when we consider non- 
perturbed and perturbed MCs (Fig. Table 0. However, our 
detailed analysis permits us to distinguish MCs that are in differ- 
ent conditions. The MCs in interaction have a broader range of ( 
(Fig.0) since they are strongly affected by the interaction, and 
they are observed at different times of the interaction process. 
The non-perturbed MCs logically have the smallest dispersion 
in £ (Fig. 0). He nce, in addition to the dispersion in £ found 
from fhe lDu etaO d2010l) results, we find that different environ- 
ments or/and evolution stages are present. 
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5. Magnetic clouds in interaction 

5.1. Evolution of MCs in a complex solar wind 

Some MCs travel in a structured SW formed by plasma coming 
from different parts of the solar corona. These plasmas have dif- 
ferent properties, for instance in terms of plasma speed and mag- 
netic field, so they are interacting. For example, a rarefaction re- 
gion is formed when a fast S W precedes a slow SW, while a com- 
pression region is formed when a fast SW overtakes a slow SW. 
Moreover, especially around the maximum of the solar cycle, 
several ICMEs could be ejected from the Sun at similar times, 
then interact later on. When a MC travels through a structured 
SW, the encountered structures affect its evolution. 

The group of perturbed MCs, analyzed above, are interpreted 
as examples of MCs affected during their travel from the Sun 
by the encountered SW structures. However, since the analyzed 
data mostly correspond to events that have interacted in the past, 
there is no direct evidence of the interaction in the in situ data. 
The previous interactions of these MCs can only be suspected 
based on an overtaking flow and a typically lower ( value. For 
some other MCs, the interaction is ongoing in the observations. 
In this case, the MC evolution depends on the type of interac- 
tion (e.g. on the external structure involved, the relative velocity, 
the magnitude of its plasma parameters, and the field strength) 
as well as the time elapsed since the beginning of the interac- 
tion. This implies that the study of MCs in interaction needs to 
be done on typically a case by case basis. We include the MCs 
that are strongly interacting at the time of their in situ observa- 
tion, in a third group (called simply interacting). In the next two 
subsections, we analyze two sub-groups (called Ib and Imc) of 
these interacting MCs. 




50 100 



500 
• 450 



?^V-..',-*l*V. 




50 100 150 200 250 300 

Time (hours), after 7 Jun 1 999 

Fig. 9. Example of two MCs in interaction. Magnetic cloud 44 is 
a large MC with a strong magnetic field; it is interacting with the 
smaller MC 26 (the vertical dashed lines define the MC bound- 
aries). The magnetic field rotation is mostly in the longitude of 
the field (O b ) for MC 44 and mostly in the latitude of the field 
(9b) for MC 26 wich implies that the MC axis are almost orthog- 
onal to each other. Nevertheless there is only a weak change in 
the magnetic field between the outward branch of MC 44 and the 
inward branch of MC 26. The MCs are separated by a sheath-like 
region with a magnetic field of different orientation. The bottom 
panel shows that MC 26 is not expanding. A linear least squares 
fit of the velocity is indicated by a green line in the time interval 
where an almost linear trend is present in each MC. 



5.2. MCs with stronger B-field nearby 

In the first sub-group, Ib, MCs have a strong magnetic field in 
their surroundings, that is typically stronger by a factor of above 
1.5 than the mean field strength inside the MC. However, this 
sub-group contains only cases where the nearby magnetic struc- 
ture does not have the characteristics of a MC. 

Magnetic clouds 2 1 and 46 both have a strong magnetic field 
within a large front sheath, which is larger than the MC size. 
Both have a low f value (0.26 and 0.39, respectively). The MCs 4 
and 27 both have a strong magnetic field behind them within an 
extended region, which is larger than the MC size. Both MCs 
have field strengths that increases from the front to the rear, par- 
ticularly MC 27 since its B-profile is unusually increases in an 
almost linear way from the front to the rear by a factor of about 
two. MC 4 has a relatively small £ value (0.54), while MC 27 
is in compression (£ = -0.81). The MCs 10 and 33 also have 
a strong magnetic field behind them, while at the difference of 
both previous MCs, their internal field strength has several rather 
unusual structures, which indicates that they are strongly per- 
turbed. They both expand less than average as £ = 0.36 and 
0.33, respectively. 

We conclude that the above six MCs are compressed by the 
surrounding strong magnetic pressure so that they expand less 
rapidly than the average rate (even one MC is in a strong com- 
pression stage). However, there is one exception, MC 1 that also 
has a strong magnetic field behind (in a region more extended 
than the MC size) but it is in over-expansion as £ = 1.45. This 
MC is likely to be in a late interaction stage where the inter- 
nal pressure, build up previously in the under-expansion stage, 



becomes too high and is able to drive an over-expansion (see 
Sect.|43]i. 

5.3. Interacting MCs 

When, after the solar eruption of a magnetic cloud, another one 
is ejected having a more rapid velocity and in a similar direction 
to the first one, an interaction between the two MCs is expected 
near their encounter time. Depending of the relative orientations 
of the magnetic fields of the two interacting MCs, magnetic re- 
connection can develop during the interaction. Previous stud- 
ies of these processes have been performed using observations 
of individual studied cases f or both favorable c onditions to the 
magnetic reconnection (e.g., Wang et al. 2005a), as well as non- 
reconnection conditions (e.g.. iDasso et al.l l2009). The dynami- 
cal evolution during the interaction of two magnetic clouds have 
been also studied numerically using MHD simulations (e.g., 
iLugaz et "ail 120051 : IXiong et al.ll2007h . These studies have typi- 
cally found that a linear velocity profile is reached in the late 
evolution phase when the two MCs are traveling together. 

We analyze these interacting MCs, with at least another MC 
in their vicinity, in the second subset of interacting MCs (Imc)- 
An example is shown in Fig. [9] The large MC 44 has a typical 
expanding rate (f = 0.95), while the smaller MC 26 located be- 
hind has no significant expansion (£ » 0.). A second couple of 
MCs have similar magnetic-field strengths (< B 1 nT), and 
the front MC 18 is in slight over-expansion (£ = 1.24), while 
behind, MC 19 is in under-expansion (£ = 0.45). A third cou- 
ple of MCs is formed by MCs 15 and 42. The MC 15 was also 
detected by ACE spacecraft at 1 AU. Its expansion rate was clas- 
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sical (f = 0.74. lNakwacki et alJl201 ll) . However, when the same 
MC is measured at 5.4 AU, it is strongly overtaken. In particular, 
its outward branch is compressed (< B >^ 1.2 nT compared to 
< B >w 0.5 nT in its inward branch) and in under-expansion 
(£ ~ 0.57 compared to ( — 0.67 in its inward branch). This 
MC 15 is overtaken by MC 42, which has a much stronger field 
(< B >k 3 nT) and is in strong over-expansion (£ = 2.35). 

Finally, the data set contains three MCs traveling close by. 
The three MCs (7,41,8) have similar magnetic field strengths, 
and this triad was identified as a single one ICME by 
iFunsten et ail d!999h with an interval of counter-streaming elec- 
trons covering the full duration of the triad. All MCs are in over- 
expansion = 3.2, 1.94, 1.5, respectively), while the triad taken 
as a whole has an expansion rate closer to the mean value of all 
the MCs ss 0.8). Since this last £ value is very different than 
the f values found within each MC, this triad of MCs is far from 
a relaxation state with a linear profile for the full time period 
(in contrast to the results of numerical simulations). At the same 
time, the three couples of MCs analyzed above also have signif- 
icantly different £ values as illustrated for one couple in Fig. [9] 

The above three pairs and one triad of interacting MCs rep- 
resent a variety of different cases. Their contrast with the simple 
evolution found in MHD simulations among similar flux ropes, 
is as follows. As with an overtaking fast flow, the MC in front 
is first compressed, so it s expansion rate decreases (Luga z~et alJ 
120051: IXiong et alj|2007l) . The MCs 18 and 19, which have simi- 
lar field strengths and sizes, are the closest to these simulations. 
However, the f ollowing MC is in under-expansion w hile it is the 
leading MC in dLugaz et al.ll2005HXiong et alj2007l) MHD sim- 
ulations. 

When the magnetic reconnection is inhibited between the 
flux ropes (by nearly parallel magnetic fields), later on in the 
MHD simulations the couple of MCs travel together forming 
only one expanding structure, i.e. wit h a nearly linear v elocity 
profile across both MCs (see Fig. 4 of IXiong etaT1l2007l) . Such 
a case was not found in our data set. 

Other MHD simulations have shown a more complex veloc- 
ity pattern, for example in the case of th e interaction of three 
flux ropes (see Fig. 6 of lLugaz et al.ll200"7l) . We still do not have 
a complete view of the interaction of MCs of various sizes, ori- 
entations, and magnetic helicities to help us interpret our obser- 
vations of interacting MCs. We also have no information about 
the history of the interaction, which is crucial to improve our 
understanding of the local measurements since the above MHD 
simulations have shown that the interaction is a time-dependent 
process (even in the simplest interaction cases). Part of the his- 
tory can be provided by other data sets such as those provided 
by heliospheric imagers and interplanetary radio data for a more 
recent time period (with the STEREO spacecraft). 

6. Summary and conclusions 

Our main goal in this paper has been to investigate how the MC 
properties, in particular their expansion rates, evolve in the outer 
heliosphere based on previous studies of the inner heliosphere 
and at 1 AU. During the travel from the Sun to the location of 
the in situ observations, the MC field may become partially re- 
connected with the overtaking magnetic field, as indicated by the 
presence of a back region in a large fraction of MCs. We have 
defined the MC boundaries to retain only the remaining flux rope 
part in order to analyze only the region that is not mixed with the 
SW. 

Each observed MC has its own properties, with different 
temporal profiles of its magnetic field and its plasma parame- 



ters. Nevertheless we have found similarities between MCs that 
have allowed us to define three MC groups. The MCs within 
the first group have a radial velocity profile that ia almost lin- 
ear with time, which implies that there is a self-similar expan- 
sion as theoretically expected, hence this group is referred to as 
non-perturbed (see examples in Fig. [TJ- In the second group, the 
temporal velocity profile can deviate significantly from linearity, 
hence they are perturbed (see examples in Fig. |2}. Finally, we 
considered separately the MCs that have an extensive magnetic- 
field structure nearby, which are defined as the group in interac- 
tion (one example is shown in Fig. [9]). 

The perturbed MCs have an average size that is smaller than 
those of the non-perturbed MCs (Fig. |4^), while the opposite 
trend is true for the magnetic field strength, proton density, and 
temperature (Fig. |4j),cj5^). All of these results are consistent 
with perturbed MCs expanding less with solar distance than non- 
perturbed MCs, a property that has indeed been found indepen- 
dently from analyses of the in situ proton velocity (Fig. [6j. For 
the proton B, we found no significant difference between these 
two MC groups (Fig.|5j)). 

The expansion in the radial direction (away from the Sun) 
is characterized by the non-dimensional parameter defined by 
Eq. (O from the in situ velocity profile. When the observed ve- 
locity profile is linear with time, the MC size has a power-law 
expansion such as S ~ Tfi [Eq. © and following text]. For non- 
perturbed MCs, £ has a narrow distribution = 1.05 + 0.34). 
This is slightly larger than the values found in both the inner 
heliosphere (0.91 ± 0.23) and at 1 AU (0.80 ± 0.18). It is also 
larger than the global expansion exponent found from the analy- 
sis of the MC size versus helio-distance (0.79 + 0.46). This lower 
exponent for the size could be due to the progressive pealing 
of the flux rope by reconnection or/and to a temporal evolution 
(Sect. SHI. 

The expansion rate of perturbed MCs, £ = 0.28 + 0.52, is 
on average significantly lower than for non-perturbed MCs, in 
quantitative agreement with previous results in the inner helio- 
sphere. The rate £ is also more variable within the group of per- 
turbed MCs, such as in the inner heliosphere, with even one MC 
in over-expansion. These results agree with a temporal evolution 
of the expansion rate during the interaction with an overtaking 
flow as found in MHD simulations. 

The MCs in interaction with a stronger magnetic field region 
have the largest variety of expansion rates f = 1 .00 + 0.93. Some 
of these MCs are in a rapid expansion stage (up to a factor of 
three bigger than for non-perturbed MCs), while one MC is even 
in compression. This wide variety of expansion rates is partially 
present in MHD simulations of interacting flux ropes or a flux 
rope overtaken by a fast stream: there is first a compression at the 
beginning of the interaction, followed by an increasingly high 
expansion rate with time that could become an over-expansion 
much latter on (with a magnitude that remains to be quantified). 
In all cases the properties of MCs that area in the process of 
interaction are difficult to interpret both because the past history 
is unknown and because they have a larger variety of properties 
than in MHD simulations. 

Appendix A: List of studied MCs 

The studied MCs are listed in Table IA.1I We added the MCs 
41 to 46, which were not present in the list of Rodriguez et al.l 
(120041) . 
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Table A.l. List of MCs and their main properties. 
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0.83 


0.10 


0.14 


60 


1.45 


2 


15-Nov-1992 23:18 


perturbed 


-20 


5.2 
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0.91 


0.04 


0.12 


82 


0.50 
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non-perturbed 


-32 


4.6 
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0.98 


0.05 


0.14 


79 
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interacting 


-52 


3.6 
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0.51 
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0.17 


0.24 


74 


0.54 
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perturbed 


-22 


1.4 
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55 
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non-perturbed 


24 


4.6 
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0.11 


67 


0.89 
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interacting 


21 


4.6 
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0.91 


0.10 


0.34 
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interacting 


20 


4.7 


512 


0.29 


0.94 


0.14 


0.31 


65 


1.50 
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perturbed 
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4.7 
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0.36 
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40 


0.41 
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interacting 
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439 
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0.19 


0.21 


80 


0.36 
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17-Aug-1997 09:59 


perturbed 


6 


5.2 


361 


0.73 


1.47 


0.25 


0.08 


79 


0.37 


12 


30-Aug-1997 19:26 


perturbed 


5 


5.2 


392 


0.35 


1.58 


0.40 


0.18 


71 


1.26 


13 
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non-perturbed 


2 


5.3 


387 


0.48 


1.04 


0.29 


0.25 


80 


1.29 


14 
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perturbed 
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0.05 
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interacting 
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89 


0.99 


24 
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0.17 
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40 


18-Jul-2002 06:51 
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41 
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55 
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42 
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interacting 
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a Number identifying MCs. 

b t c is the time of closest approach from the MC axis. 

The MCs are separated in three groups: non-perturbed, perturbed, and in interaction (denoted interacting). 

d 6 and D are the latitude and solar distance. 

c V c is the velocity at the closest distance from the MC axis and S is the MC size. Both are computed in the radial direction away from the Sun 
(R). 

f < B >, < N p > and < B p > are the average over the flux rope of the field strength, the proton density and the proton B, respectively. 

E y is the acute angle between the MC axis and the radial direction (R). 

h ( is the unidimensional expansion rate [Eq. y}]. f < means a MC observed in compression stage. 
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